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Abstract. We investigate the dark halo lens mass function (MF) for a wide class of spheroidal non singular 
isothermal models comparing observed and observable microlensing quantities for MACHO observations towards 
LMC and taking into account the detection efficiency. We evaluate the microlensing observable quantities, i.e. 
observable optical depth, number of events and mean duration, for models with homogenous power -law MF 
changing the upper and lower mass limits and the flattening of the dark halo. By applying the simple technique 
of the inverse problem method we are then able to get some interesting constraints on the slope a of the MF and 
on the dark halo mass fraction / made out by MACHOs consistently with previous results. 
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1. Introduction 

The usefulness of gravitational microlensing as a 
tool to investigate the structure of the galactic dark 
halo is now firmly established. Since Paczynski's 
seminal paper (Paczynski 1986), several collabora- 
tions have searched for microlensing events towards 
LMC ([Renault et al. 1997|; |Alcock et al 2000a| ), SMC 
( lAlcock et al. 1997c[ [Afonso et al. 1999| ) and M31 
( [Ansari et al. 1997| ). While the presence of MACHOs 



with radius equal to the Einstein radius, defined as 



(Mollerach fc Roulct 1997) 



Re 
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(Massive Astrophysical Compact Halo Objects) has been 
witnessed by the microlensing events already found, many 
questions about their nature and the structure of the 
dark halo are still open. 

One of the main quantities that characterise the mi- 
crolensing properties of a particular MACHO popula- 
tion is the optical depth r. This quantity is simply the 
number of lenses inside the microlensing tube, a cylin- 
der whose axis is the line of sight to the source and 



where m is the mass of the lens ^] (with p, = tu/Mq), D s 
(Dj) the distance to the source (the lens), s = Dd/D s and 
we have posed r E = \J AGMqD s / 'c 2 . In its simplest form , 
t is defined as ( Paczynski 1986| ; Mollerach fc Roulet 1997 ; 
Jetzer 1998|) 



^-p{s)ds 
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p(s)s(l - s)ds , (1) 



having supposed (as usually) that the halo extends till 
the source |, while p(s) is the MACHOs mass density. 
Optical depth is obtained from observations using the for- 
mula ( |Alcock et al. 1997b| ) 
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In this paper we use the terms lens and MACHO as syn- 
onymous, even if the lens may be also a visible star. 

2 This is a reasonable hypothesis if the source is in Magellanic 
Clouds, but it is not if the target of observations is M31. 
However, in this paper we limit ourselves to observations to- 
wards LMC only. 
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being E the total exposure in star - years (which is simply 
the number of monitored object multiplied by the total pe- 
riod of observations), ts = ^Re/v±_ the Einstein diameter 
crossing time (with v± the MACHO' s velocity transverse 
to the line of sight), s(tE,i) the detection efficiency for the 
i-th event and the sum is over the observed events. 

Having estimated T meas from observations with Eq. (|J) 
and calculated the predicted optical depth for an assumed 
model of full MACHO halo from Eq.©, an easy way to 
evaluate the halo mass fraction composed by MACHOs is 
to compute the ratio between these two quantities to get 
/ = T meas /r. However, Alcock et al. (1997a) have pointed 
out that Eq. (Q) is not a measure of the total optical depth, 
but only of the optical depth of events which fall in the 
range of timescales for which > 0. Following this 

remark, Kerins (1998) has introduced the concept of ob- 
servable optical depth defining it as 



Toble t 



e(* E )>0 



, , dY 
£(tE)tE-^—dtE , 

aiE 
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where dY/dtE is the differential rate for the assumed 
halo mode l, being the rate Y the number of eve nts per 
unit time ( Mollerach fc Roulet 1997 ; Jetzer li 



Then, 

the halo mass fraction / is more correctly estimated as 
Tobs d/joblej being Tobsd the directly observed optical depth, 



i.e. (iKerms 1998[ ) 



Tobsd ~ 4E 2_^ 
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Using these quantities, Kerins has shown that the com- 
parison between observed and obervable optical depths 
may give only a lower limit on /, while no upper limit 
may be obtained without a knowledge of the full distri- 
bution function of the halo model and of the MACHO 's 
mass function (hereafter MF), i.e. the number density of 
MACHOs with mass in the range (/i, \i + d/i). The re- 
sult obtained by Kerins was proved under the hypothesis 
that all the lenses have the same mass and assuming a 
standard cored isothermal sphere for the dark halo mass 
density. In this paper we generalize the calculation of T Ue 
and analyze a general class of homogenous MF. Several 
studies have been made by many authors to determine 
the MF function by microlensing data, but they are essen- 
tially devoted to the MF of bulge lenses ([Zhao et al. 1995 ; 
Han & Gould 1996; Grenacher et al. 1999[ ) and the work 
is still in progress. Mao & Paczynski (1996) considered 
simplified toy models and a power - law MF and estimated 
that a reliable determination could be achieved only if we 
had 100 or more events. Their results, obtained under the 
assumption that MACHO 's spatial distribution and kine- 
matics were known, have been confirmed by Markovic & 
Sommer - Larsen (1997) who have also studied the effect 
of changing halo model on the estimated average mass of 
lensing objects. All these studies are based on statistical 
methods, that is why there is need of a large number of 
events to reduce the error on parameters determination. 
In this paper we use a different technique to get useful 



informations on the lensing objects MF. Since microlens- 
ing observable quantities, i.e. number of events, observable 
optical depth and mean duration, depend on the dark halo 
model and on the MF, we may express these quantities as 
functions of some models parameters and of the slope of 
the MF itself. Then a comparison between theoretical ex- 
pectations and observed quantities will help us to recover 
the values of the model's parameters simply imposing that 
theory and observations are in good agreement. This is 
what we call the inverse problem in microlensing. To be 
meaningful our analysis must take into account also the 
detection efficiency and that is why our theoretical expec- 
tations have to be corrected for this effect too. In the fol- 
lowing we will show how this is possible and will get some 
interesting constraints on the dark halo mass fraction / 
composed by MACHOs and the slope a of the assumed 
homogenous power - law MF for a wide class of spheroidal 
models. By the way our method is not able to escape the 
problems connected to the low number of observed events 
and actually our results are strongly affected by errors 
which do not allow us to constrain the slope a of the MF 
in a narrow range. Were the errors be reduced by increas- 
ing the number of observed events, our method should be 
able to narrow the uncertainties on the MACHOs MF. 

There is one possible source of systematic error con- 
nected to our analysis. In the previous discussion, we have 
implicitly assumed that all the observed events are due to 
MACHOs in the Milky Way dark halo, but it should be 
considered also the possibility that at least part of the 
events are due to LMC self lensing. This hypothesis has 
been suggested in many papers, but the recent analysis by 
the MACHO group of the spatial distribution of the events 
across the observed LMC fields has shown that this latter 
is not compatible with the proposed models of self lensing. 
As a further test, MACHO has also analyzed the CMD 
(Color Magnitude Diagram) of the sources of their first 
eight events to see if they reside in the LMC disk or behind 
it. Even if the sample is too small to get a definitive result, 
the hypothesis that all th e eight events arc d ue to halo 
lenses is slightly favoured ( Alcock et al. 2000b ). However, 
there is still the possibility that there are no MACHOs at 
all in the dark halo and that the observed events are due 
to the LMC dark halo. In fact, the maximum likelihood 
analysis of the MACHO collaboration has shown that this 
hypothesis is not completely excluded, even if the needed 
LMC model should be somewhat extreme. Besides, a re- 
cent paper by Alves & Nelson (2000) has shown that the 
LMC rotation curve and the data on the kinematics of the 
carbon stars in LMC are best fitted by a model composed 
by a flared and warped disk and no dark halo. If this result 
will be confirmed by future works, it will strenghten our 
assumption that the observed microlensing events are due 
to MACHOs in Milky Way dark halo. In the meantime, 
we are constrained to consider the effects that a possible 
contamination by self lensing should have on our results. 
We will discuss how this will affect our results doing some 
qualitative tests. 
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In Sect. 2 we introduce and evaluate the microlensing 
observables, i.e. the number of events, the observable opti- 
cal depth and the mean duration, taking into account the 
detection efficiency for models with isotropic maxwellian 
transverse velocity distribution and homogenous power - 
law MF and assuming a wide class of spheroidal non sin- 
gular isothermal models. The model parameters and the 
way we fix some of them is detailed in Sect. 3 where we 
illustrate the models we have chosen to explore. Sect. 4 is 
devoted to the analysis of the different models with the 
simple technique of the inverse problem method to get 
constraints on the slope of the lens MF and the dark halo 
mass fraction composed by MACHOs. How the self lensing 
could affect our results is discussed qualitatively in Sect. 
5, while a final discussion of the results is presented in the 
conclusions. 

2. Microlensing observables for spheroidal models 
with power -law MF 

To apply the method we have outlined in the introduction, 
we have first to calculate the observable quantities in mi- 
crolensing to which we will compare the observed ones 
towards LMC. These quantities are the number of events, 
the observable optical depth and the mean duration. In all 
this calculation, we will take into account the detection ef- 
ficiency in order to be sure that the comparison between 
predicted and observed quantities is meaningful. 

The starti ng point is the general expression of the dif- 
ferential rate ( Dc Rujula et al. 1991 ) 



dT = 2D s r E v ± f(v ± )[p,s(l 



s)\ 1 —dpdv±ds 



dp 



(5) 



In Eq. (JSJ) , f(v±) is the transverse velocity distribution and 
dn/dp(s, p) is the halo lens MF; in the reasonable hypoth- 
esis of homogeneity, this latter may be factorized as : 



dn dn p(s) dn 

H{s) x— („) = — x— 0*) 



(6) 



being p Q the local mass density of the model and dno I dp 
the local MF. Changing variable from v± to t E and inte- 
grating we get 



dT 

dt E 



SD s r 3 E 

1 E 



>(1 - s)] 3/2 ds 



/■'-' 



3/2 dn Q ^ 
dp i 



(7) 



being {pi,p u ) the lower and upper limit for the mass of 
MACHOs. From Eq.(0) one immediately see that to go 
on further we need to assign the transverse velocity dis- 
tribution, the local MF and the mass density of the halo 
model. 

Since we do not consider anisotropy in the velocity 
space, we may assume the following maxwellian distribu- 
tion of the transverse velocities 
2w ± 



/Ox) 



exp{-v 2 Jv 2 H } 



(8) 



where vr is the velocity dispersion which we fix as 210 
km/s. 

With regard to the local MF, it is usual to assume that 
all the MACHOs have the same mass which means that 
dno /dp is a (S-Dirac centred on the common mass. This 
is just a first approximation: it is worthwhile to explore 
different possibilities. As a generalization we consider the 
case of a homogenous power -law MF for the MACHOs, 
i.e. we assume 



dno 
dp 



(p) = C(a)p- a 



(9) 



being C{a) a normalization constant fixed such that 



C(a) 



this gives 



p 1 - a dp 



I 2-l' a 'A- a (a ^ 2) 



In fi u —In Hi 



(a 



(10) 



(11) 



Concerning the mass density of the dark halo we will 
restrict our analysis to a class of spheroidal non singular 
isothermal models whose density distribution is given by 



p(R,z) = p { Q } 



R + R c 



qsm- 1 eR 2 + z 2 /q 2 +R 2 c 



(12) 



t t (s) . 

where e = s/1 — q with q the halo flattening, p & is the 
local mass density^ for the spherical case, Rq is the galac- 
tocentric distance of the Sun and R c is the core radius. 

Inserting now Eqs.(|), (§) and (||) into Eq.fljj) and 
expressing (R, z) in terms of (s, I, b) (with (I, b) galactic 
angular coordinates of the target), we finally get 

d£_ _ ^^D s C{a)r 2 E {a - l \Rl + R 2 ) 
dt E 



2(a-2) Ma-\) 
U H l E 



1 [*(!-*)] 



Q-l 



G(s, t E )ds 



(13) 



Jo As 2 +Bs + C 
being 

A = D 2 (cos 2 b + q~ 2 sin 2 b) , 
B = —2D s Rocosbcosl , 
C = R 2 + R 2 , 

and we have defined for sake of shortness 

G(s, t E ) = G(3 - a, 0{s, t E )m) - G(3 - a, f3(s, t E )p u )(U) 

with 



(3(s,t E ) 



4r|s(l - s) 



vU 



, G(o,C) 



HE 



C 



t a_1 exp {-t}dt 



a Note that the pq previously introduced and Pg' are equal 
only for models with q = 1; in general, the link between these 
two quantities is given by: p Q = Pe to -i e - 
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Given p(R, z) it is quite straightforward to get all the 
observable quantities we need in the following. The first 
one is the number of observable events, which is simply 
given by 



jyoble _ E 



e(t E )>0 



£\t E )——dt E 
at E 



(15) 



Introducing the expression given by ( |l3| ) into Eq. (|l5|) one 

gets 



N°l le = K ev (a)I ev (a; model) 
having posed 



K ev (a) 



2 ( ~ 2a -^D s C(a)E: 



2(a-l) 
E 



(Rl + R c ) 



2(a-2) 
H 



(16) 



(17) 



lev {a; ma 



del) = [ 

Je(t E )>0 



e{t E )t 



-2(a-l) 



E 



dti 



(4e)>( 

1 [s(l-s)]"- 1 
o As 2 + Bs + C 



G(s,t E )ds 



(18) 



Also, from Eq.([l3|) and definition (||), we get the following 
expression for the observable optical depth: 



T ue = K T (a)I T (a; model) , 
being now 

2(^-3) 7TDsC (a)r 2 E {a - 1) (Rl + R 2 C ) 



K T (a) = 



' E 
2(u-2) 



(19) 



(20) 



I T (a; model) = / s{t E )t 

Je(t E )>Q 



3 ~ 2a dt E 



As 2 + Bs + C 



G(s, t E )ds 



(21) 



Now it is not difficult to get the predicted mean duration; 
in order to take into account the detection efficiency, we 
define this quantity as 



< t E >= 



1 



old, C 



AToble , 

Iy ev Je(t E )>0 



t E dN°° 



(22) 



This definition is nothing else but a straightforward gen- 
eralization of the usual one (Jetzer 1998) to which it re- 
duces in the case of perfect detection efficiency, i.e. when 
e(t E ) — 1 V t E . It is not surprising, then, that the follow- 
ing relation holds: 



Table 



AE 



<t E > N { 



oble 



which is immediately demonstrated. Using such a relation 
and Eqs.rtll) and (19), we finally get 



< t E >= model) . 



(23) 



In Eqs.©, (£§) and Q we have indicated the unknown 
parameters which these quantities depend on. These are 



(s) 

Table 1. Models main characteristics. The parameter p Q 
is in units of 10 -3 M & pc~ 3 ; M50 is the total mass inside 
50 kpc in units of 10 11 Mq and v c (oo) the asymptotic 
rotation velocity in km/s. 



q 




M 50 


Vc(oo) 


0.3 


6.7 


3.90 


197 


0.5 


7.2 


4.02 


201 


0.8 


8.1 


4.21 


208 


1.0 


8.7 


4.32 


212 



essentially the slope a of the MF and the parameters 
needed to specify the characteristics of the mass model. 
We have decided to not consider the mass limits as un- 
known parameters, although they are not very well con- 
strained. In fact, in the following analysis we retain as 
fixed these two quantities and investigate the dependence 
on the remaining unknowns. 

3. Exploring the models parameters space 

The integrals in Eqs.© and © may be evaluated only 
numerically, but to do this we need to fix all the parame- 
ters entering the mass density in Eq.(p"2|) (i.e. q, p^', R a , 
i? c ), the lower and upper mass limits (/i;, p u ) and also the 
slope a of the MF which is just the parameter we want to 
determine. 

We have fixed Rq = 8.0 kpc and R c — 5.6 kpc (as 
usual in literature). Next we have to choose the halo flat- 
tening q, but the constrains on it are really very poor. 
Different kinds of analysis with different techniques give 
very different results (for a review see, e.g., Rix 1996 and 
[Sackctt 1999 ) and there is no general agreement on which 
value is the best one to use. Since there is such a large un- 
certainty, we have decided to consider four values of the 
halo flattening, repeating our analyses for models with 
q = 0.3,0.5,0.8, 1.0 in order to test if this parameter has 
some effect on the results. The last parameter to fix is 
the local mass density for the spherical case Pq\ This 
is fixed such that the local rotation velocity is equal to 
the observed value of v c (Rq) — 200 km/s. To this aim we 
have also considered the contributions of the bulge and 
disk. We modelled these components as in Mera et al. 
(1998): the bulge is treated as point like with total mass 
1.23 x 10 10 Mq and the disk as a double exponential with 
scale length Rd = 3.5 kpc and local surface mass density 
S Q = 52 M & pc~ 2 . In Table 1 we give the values of the 
models parameters together with the total mass inside 50 
kpc and the asymptotic rotation velocity. We would like 
to note that the predicted values of these latter quantities 
are in good agreement within the errors with the recent 
measurements (see, e.g., Wilkinson fc Evans 1999] ) for all 
the considered models. 

In order to obtain the observable in term of a we have 
numerically integrated Eqs.(18) and (^l|) for many values 
of a and then interpolated the results to get I ev and I T as 
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Table 2. Values of the coefficients entering the numerical approximations given in Eqs.(|25|) and ( pq ) for models Al 
(first half) and A2 (second half). 



q 






C-ev 






6r 




d T 


0.3 


0.08 


0.22 


-8.83 


3.42 


0.13 


-0.09 


-8.67 


7.00 


0.5 


0.08 


0.23 


-8.84 


3.94 


0.13 


-0.09 


-8.67 


7.57 


0.8 


0.08 


0.24 


-8.84 


4.26 


0.13 


-0.09 


-8.67 


7.92 


1.0 


0.08 


0.24 


-8.84 


4.36 


0.13 


-0.09 


-8.67 


8.02 


0.3 


0.17 


0.35 


-11.33 


7.48 


0.29 


-0.47 


-10.67 


11.92 


0.5 


0.17 


0.38 


-11.35 


7.98 


0.29 


-0.47 


-10.67 


12.47 


0.8 


0.17 


0.39 


-11.36 


8.29 


0.29 


-0.46 


-10.68 


12.81 


1.0 


0.16 


0.40 


-11.37 


8.38 


0.29 


-0.46 


-10.68 


12.91 



Table 3. The same as in Table 2 but for models Bl (first half) and B2 (second half). 



q 


(lev 


bev 


Cev 




a T 


b T 


Ct 




0.3 


0.008 


0.16 


-8.75 


3.38 


0.02 


0.10 


-8.81 


7.02 


0.5 


0.008 


0.16 


-8.75 


3.89 


0.02 


0.10 


-8.81 


7.58 


0.8 


0.009 


0.16 


-8.74 


4.21 


0.02 


0.10 


-8.81 


7.93 


1.0 


0.008 


0.16 


-8.75 


4.31 


0.02 


0.10 


-8.81 


8.04 


0.3 


0.05 


0.42 


-11.23 


7.44 


0.12 


-0.01 


-11.02 


11.98 


0.5 


0.05 


0.43 


-11.23 


7.93 


0.12 


0.004 


-11.02 


12.53 


0.8 


0.05 


0.43 


-11.24 


8.24 


0.12 


0.007 


-11.02 


12.87 


1.0 


0.05 


0.43 


-11.24 


8.34 


0.12 


0.008 


-11.02 


12.97 



functions of a itself. We have then repeated this procedure 
changing the values for the halo flattening q and the mass 
limits in order to investigate a wide class of halo models, 
each one labelled with a code given as follows. We named 
Al, A2 models with \n = 0.001 and \x u = 0.1 and fi u = 1.0 
respectively, and Bl, B2 models with /i; = 0.01 and fi u = 
0.1 and /j, u = 1.0 respectively. Then we add a letter to 
indicate the halo flattening with the following conventions: 
a -> q = 0.3, b -> q = 0.5, c -> q = 0.8, d -> q = 1.0. 
So, e.g., the model labelled A2c has: m = 0.001, \i u = 1.0, 
q = 0.8. Thus, we consider sixteen different models in the 
same class. 



Before going on, we would like to discuss how 
we have chosen the mass limits (ni,/i u ). Concerning 
the upper limit, it has been known for a long time 
now ( |Gilmore fc Hewctt 1983 ) that hydrogen - burning 
stars cannot provide the majority of the halo dark 
matter. Numerous recent studies flBahcall ct al. 1994 ; 
Hu et al. 1991 ; |Graff fc Freese 1996] ) put an upper limit 
of at most 4% on the dark halo density contribution of 
hydrogen - burning stars. For an old, metal- poor popu- 
lation this means that stars with mass between 0.1 M Q 
and 0.8 Mq give no significant contribution to the dark 
matter in the galactic halo. These evidences suggest to fix 
fx u = 0.1, but we have decided to consider also models 
with fi u = 1.0 for the reasons we are going to explain. 
On one hand, MAC HO QAlcock et al. 2000aD and EROS 
( Renault et al. 1997 ) results indicate that the most likely 
MACHO' s mass is 0.5 Mq. On the other hand, Kerins 
(1997) has shown that MACHOs may reside in a popu- 
lation of dim halo globular clusters comprising mostly or 
entirely low -mass stars just above the hydrogen - burning 



limit. For the case of the standard halo model, this sce- 
nario is consistent not only with MACHO observations, 
but also with cluster dynamical constraints and number 
counts limits imposed by twenty HST fields. Further sug- 
gestions of the possible existence of MACHOs with mass 
~ 0.5 Mq come from the study of the double quasars 



variability (Koopmans & dc Bruyn 2000). All these stud- 
ies have led us to consider also models with fi u = 1.0. 

Fixing the lower limit /!/ is not an easy task too. De 
Rujula et al. (1991) have shown that a lower limit for 
the mass of MACHOs is 10~ 7 Mq, but this does not im- 
ply that objects with a mass ~ 10~ 7 Mq really exist. 
Actually, MACHO and EROS search for short duration 
events pose strong constraints on their contribution to the 



halo mass budget ( 


Alcock et al. 1996|; 


Alcock et al. 1998 


Renault et al. 1998 


). Following these works, we have cho- 



sen two values for \i\ given by 0.001 and 0.01 respectively. 
In all our analysis we are assuming that the MF is the 
same in the mass range (/i;,/z u ) Mq, i.e. that the slope a 
does not change in this range, which seems quite reason- 
able as a first approximation. 

To integrate Eqs.® and @ we need the detection 
efficiency of the MACHO collaboration for their moni- 
toring campaign towards LMC since in our analysis we 
will use their results of the first 5.7 years of observations. 
This function has been carefully e valuated by the MACHO 
group itself (Alcock ct al. 2000b), but they give no ana- 



lytical formula for it. That is why we have built up an 
approximated expression of e(t E ) interpolating the data 
taken from Fig. 5 of Alcock et al. (2000a), obtaining 



e{ti 



0.108 exp{-0.02(logt£;) 3 - 8 }(logi E ) 2 - 4 
0.41exp{-14.7(logt£; - 2.56) 2 } 



(24) 
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being log the decimal logarithm; the first expression holds 
for 2 d < t E < 300 d and the second for 300 d < t E < 



900 d. Eq.(24) differs from the measured e(t E ) less than 
10% in the range examined, the error being larger for 
events lasting more than 900 days. This is not a problem 
since the thirteen observed events which we use in our 
analysis last approximately from 34 to 102 days, so we are 
confident that no serious systematic error is induced by 
our approximation for e(t E ). For the same reason also the 
discontinuity in t E = 300 d has no effect on our analysis. 



Fig. 1. < t E > (a) for the model named A2c, i.e. with 
\H = 0.001, fi u — 1.0, q — 0.8. Similar plots are obtained 
for the other models. 

Now we have all we need to estimate the functions 
defined in Eqs.(p"8|) and (pl|). Without entering in details, 
the numerical integration and the following interpolation 
of the results have shown us that it is possible to write 



I ev {a) = exp{a etI a 3 + b ev a 2 + c ev a + d ev } 



I T (a) = exp{a r a 3 + b T 



d T } , 



(25) 
(26) 



where the values of the coefficients (a ev , b evi c ev , d ev ) and 
(a T , b T ,c Tl d T ) depend on the model considered. We report 
them in Tables 2 and 3. These approximations work quite 
well, the error being always less than 10 % and in many 
cases also less than 5 %. As an example, in Fig. 1 we plot 
the ratio I T /I ev =< t E > (a) for model A2c. 

4. The inverse problem to estimate (a, f) 

We show now how it is possible to use the simple tech- 
nique of the inverse problem to estimate the slope a of 
the lens MF and the dark halo mass fraction / composed 
by MACHOs for each one of the models labelled by the 
codes previously explained. We consider data coming from 
the first 5.7 years of observations towards LMC by the 



MACHO collaboration ( |Alcock ct al. 2000a[ ), limiting our- 
selves to the thirteen events selected according to the so - 
called selection criteria A. These are high S/N events and 
are spatially distributed in a way which is consistent with 
the hypothesis that they are due to lenses belonging to 
our halo and not to LMC self lensing; we will discuss of 
this problem later on. For this set of events we have 



N. 



obsd 



= 13 , r obsd = 3.5 x 10- 



< *e >obsd= 76.2 d 



where < t E > bsd is simply the average value of the dura- 
tion of the observed events. As a test of the correctness of 
this estimate we may note that (n/AE) < t E > bsd (with 
E = 6.12 x 10 7 star-years) is exactly equal to T bsd/N° v ' sd , 
as it has to be for the reasons we are going to explain later. 
We have then to estimate the uncertainties on the ob- 
served quantities. First, we consider the directly observed 
optical depth, and we simply use a method as similar as 
possible to the one proposed by the MACHO group it- 
self for a very conservative estimate of the error on T meas 
(Alcock et al. 1997a): we divide the observed events ac- 
cording to their duration t E in bin of 10 days; in such a 
way T bsd is more or less the same for events in the same 
bin and the errors are approximately poissonian. For each 
bin we estimate Ni ow = N°^ sd (bin) 
(bin) 



y/N°^ sd (bin) and 



N — N obsd 



y/N^ sd (bin) and define r 



mm and 



T max ag minimum and maximum value of t\ for the 
events in that bin (being n — ixt E /AE) . Then we estimate 



obsd 



E 

bin 



AT rr-rn 
NlowT i 



obsd 



5> 



up ' 1 



bin 



At the end we find 

Tobsd = (3.5 ±2.5) x 10" 8 . 

The error on the optical depth turns out to be so large (~ 
70%) because of the limited number of events. Let us turn 
now to the error on the number of observed events. This 
simply comes from the low statistics and may be assumed 
to be poissonian, i.e. SN^ sd /N^ sd = l/^N^ sd ~ 28%. 
Finally, the error on < t E > is obtained by propagating 
the error on N°^ sd which gives us 



8 < t E > , 



5N° : 



bsd 



< t E >obsd 



ATobsd 
ev 



The observed quantities may be now compared to the 
theoretical ones evaluated in the previous section, which 
also take into account the detection efficiency in order to 
make the comparison meaningful. To do this we must first 
remember that the theoretical quantities have been calcu- 
lated under the hypothesis that the dark halo is totally 
made out by MACHOs, i.e. with / = 1. Actually, the ex- 
act value of / is not well known: from the more recent 
observational constraints it is quite unlikely that f — 1. 
However to take into account of / is not very difficult: we 
have simply to multiply by / the expression of the differ- 
ential rate dT/dt E and consequently the ones obtained for 
N°^ le and T u e - Note that < t E > is independent on the 
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Fig. 2. In the upper panel we plot / x N°% as a function of a and /; in the lower one, we determine the region 1Z 
of the parameter space (a, /) consistent with microlensing observations towards LMC: the upper line in the second 
panel is the level curve / x N° ble — 16.61, the lower one for / x N° ble = 9.39, whilst the two dashed lines have been 
drawn intersecting the level curves with the vertical lines a — 0.79 and a = 1.50. The plot is for model A2c; similar 
plots are obtained for the other models. 



value of /. Then we have the following relations between 
observable and observed quantities: 



iwoble 

ev 

Toble = 



-- N° b v sd /f 

Tobsd/ f 



(27) 



Dividing these two equations and using the relation 
Table/ N°* le = (?r/4£) < t E >, one gets: T obsd /N° b v sd = 
(ir/AE) < t E >obsd, which may be used to test the cor- 
rectness of our previous estimate of < t E >obsd- For each 
model, we may solve the system ( |27| ) in the unknowns a 
and / and estimate them toghether with the errors con- 
nected to our analysis. Solving Eqs.(|27]) is not really useful 
since the very high error on T b s( i leads us to get no con- 
straints at all on the parameters (a, /). We may then use 
a third relation we have at our disposal, given by 



<t E > (a) =< t E >obsd 



(28) 



which is quite easy to solve numerically to get an estimate 
of the slope a of the MF. Solving Eq.(|2q), taking into 
account also the uncertainties, will give us in general more 



than one solution compatible with the microlensing data. 
A further selection can be done imposing that a must be 
in the range (0.0, 5.0) since values outside this range are 
reasonably quite unlikely^]. In this way, for each model, we 
have estimated a range (ai,a u ) for the slope of the MF 
simply requiring that 

(1 - 5t E /t E ) < t E >obsd<< t E > (a) < 



(1 + St E /t E ) < t E > obsd 
54.86 d<<t E > (a) < 97.54 d 



(29) 



being 8t E /t E the fractional error on < t E > Q bsd which we 
have previously estimated to be of order of 28%. These 
values are summarised in Table 4, where we report also 
the value of cxq corresponding to the average value < t E > 
(a ) = 76.20d. 

4 Negative values of a means that the MF is increasing which 
has never been observed, whilst values greater than 5.0 are 
quite unusual. 
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Table 4. Estimates of the slope a and the dark halo mass fraction / composed by MACHOs for different models 
labelled as explained before. M^ r is the mass in MACHOs inside 50 kpc (with M^ r ^ the value corresponding to cvq 



'50 

ao 



and /o) measured in units of 10 M@. Note that the value of ai for model B2a differs considerably from the others 
since the maximum value of < t E > (a) is only 96.2 d, which is lower than the upper limit on the < t E >obsd given 
in Eq.Q. 



Code 


{ai,a u ) 


[fmin-) fmax^j 




(q ,/ , A/ 5 b o, r (0 )) 


A2a 
A2b 
A2c 
A2d 


(0.50, 1.48) 
(0.73, 1.50) 
(0.80, 1.50) 
(0.84, 1.56) 


(0.12, 0.33) 
(0.10, 0.30) 
(0.11, 0.32) 
(0.12, 0.35) 


4.68 - 
4.02 - 
4.63 - 
5.18 - 


- 12.87 

- 12.06 

- 13.47 

- 15.12 


(1.10, 0.21, 8.19) 
(1.15, 0.19, 7.63) 
(1.17, 0.20, 8.42) 
(1.21, 0.22, 9.50) 


B2a 
B2b 
B2c 
B2d 


(0.27, 1.69) 
(0.63, 1.80) 
(0.77, 1.89) 
(0.77, 1.90) 


(0.11, 0.33) 
(0.10, 0.30) 
(0.10, 0.32) 
(0.11, 0.34) 


4.29 - 
4.02 - 
4.21 - 
4.75 - 


- 12.87 

- 12.06 

- 13.47 

- 14.79 


(1.12, 0.24, 9.36) 
(1.25, 0.18, 7.24) 
(1.34, 0.19, 8.00) 
(1.34, 0.21, 9.07) 



Having estimated a and being N°° sd known, we may 
now get a constraint also on /: on the plot fN°^ le as 
a function of a and / we make the contour levels for 
fN°* le = N°* sd - 28 % = 9.39 and fN°* le = N°* sd + 
28 % = 16.61. On this graph one has to add also the ver- 
tical lines corresponding to a = ai and a = a u (see Fig. 
2). The region 1Z of the parameter space (a, /) delimited 
by the two level curves and these two vertical lines is that 
which is consistent with the constraints on a and the ones 
coming from microlensing observations towards LMC, i.e. 
the number of observed events and their mean duration. 
From Fig. 2 one sees that it is possible to define f m in and 
fmax as the minimum and maximum value of / in the 
region 1Z and a value /o such that 

| f N^ le (a ) = 13 

[ < t E > (a ) = 76.20 d 

One has then that for each value of / G (fmin, fmax) there 
exists a value of a € (ai,a u ) such that 

( 9.39 < fN°% le (a) < 16.61 

{ . (30) 

[ 54.86 d<<t E > (a) < 97.54 d 

The values of [fmin, fmax) and /o are summarised in 
Table 4 for the different models considered^]. We also give 
in the same table the range for the mass in MACHOs in- 
side 50 kpc (indicated as M^ r ) which is easily estimated 
as / x M50 (with M50 given in Table 1) since it is this 
quantity which is most strongly constrained by microlens- 
ing observations. 

5. The self lensing problem 

We discuss now the problem of the LMC self lensing and 
show how our results depend quite weakly on the sys- 
tematics connected with it LMC. As discussed in the in- 
troduction, our analysis implicitly assumed that all the 

5 A complete list of the solutions for each model, comprising 
also the rejected ones, and further details on the numerical 
procedures are available on request to the authors. 



observed events are due to lenses belonging to our dark 
halo. Even if there are different evidencies against a dom- 
inant contribution of self lensing, it is still possible that 
two or three of the 13 events considered are due to LMC 
self lensing. This means that we have to repeat our analy- 
sis excluding from the sample those events which are not 
caused by MACHOs. Unfortunately, it is not really pos- 
sible to establish precisely what are these events, due to 
the well known degeneracy in the lens parameters which 
do not allow to determine its distance with respect to the 
observer. However, we may qualitatevely correct our re- 
sults excluding form the sample some events chosen ac- 
cording to which are more likely to be due to self lensing. 
To this aim we have used the results on the self lensing 
events timescale distribution obtained by Gyuk, Dalai and 
Griest (2000). Their analysis shows that the most likely 
duration of self lensing events is ~ 100 d; so we repeated 
our analysis for different choices of the excluded events. 
As an example we discuss here the results obtained for 
the model A2c (m = 0.001, [i u — 1.0, q — 0.8); similar 
results are obtained for the other models. As a first test, 
we have excluded from the sample the events labelled 7, 
13 and 14 in Table 7 of Alcock et al. (2000a) which are 
the longest ones. We have then < t E > b s d— 68.78 d and 
N°^ sd = 10; the uncertainties on these quantities increase 
to ~ 32% both on< t E > obsd and on N°% sd . We get 

a = 0.94 1.67 ; / = 0.07 0.24 ; (a , /o) = (1-30, 0.14). 

Even if there is a little trend towards larger values of a 
and smaller values of /, these results are quite consistent 
with the ones given in Table 4. The central value for the 
slope a is slightly larger and fo slightly lower but the 
discrepancies are not significative. We have repeated the 
same analysis in other two cases: excluding events 5, 7 
and 13 and exlcuding events 13, 21, 25 respectively. The 
results we get are consistent with each other and with the 
results for model A2c in Table 4. Finally we have also 
tried to exclude four events (events number 7, 13, 14 and 
21) from the sample; the uncertainties are of course larger 
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(~ 33% both on < ts >obsd and on N°^ sd ) while the mean 
duration is 66.07 d for 9 events. We obtain 

a = 0.99 1.72 ; / = 0.06 0.21 ; (a , fo) = (1-34, 0.12). 

These results are still consistent with those obtained till 
now for model A2c. These qualitative tests are encourag- 
ing since they suggest that our analysis is not seriously 
affected by the systematics connected to the self lensing 
whose main effect seems to be to lower the statistics. 



6. Discussion of results and conclusions 

In this paper we have analysed a wide class of spheroidal 
non singular isothermal dark halo models with a homoge- 
nous power - law MF by changing the halo flattening q and 
the lower and upper MACHO mass limits (fxi,fj, u ). Using 
the simple technique of the inverse problem method we 
have obtained useful constraints on the slope a of the MF 
and the dark halo mass fraction / made out by MACHOs. 
The results obtained are summarized in Table 4, where 
we report the estimated ranges for the slope a of the lens 
MF, for the the dark halo mass fraction / made out by 
MACHOs and for the mass in MACHOs inside 50 kpc. In 
the same Table we also report the values (ceo, fo, M^ r ^) 

(the latter being the value of M^§ r obtained with / = fo) 
which may be considered as a sort of best fit for each 
model, since for these values of the parameters the model 
is consistent with microlensing constraints independently 
on the estimates of the errors. Having in mind that for 
each value of / £ (fmim I max) there exists only one value 
of a € (ai,a u ) such that the model with these values of 
(a,/) satisfies the constraints (j30|), we may draw some 
interesting considerations from Table 4. 



1. The first striking feature in Table 4 is that there are 
no models with fi u = 0.1 (models Al and Bl in our 
notation). This happens because the constraint ( p9[) 
on < tE > (a) for these models gives values of a out- 
side the range (0.0,5.0). Thus we may conclude that 
non singular isothermal spheroidal models filled with 
MACHOs of mass less than 0.1 A/q are not consistent 
with microlensing observations towards LMC. Thus 
there exist MACHOs with mass greater than the hy- 
drogen burning limit (~ 0.1 Mq); such interpretation 
is consistent with MACHO and EROS result on the 
most likely MACHO mass (~ 0.5 Mq) and will have 
important implications on the question of their na- 
ture. Actually, if MACHOs weigth more than 0.1 M Q , 
then they should be unusual baryonic objects (e.g. very 
old white dwarfs) or non baryonic at all flSazhin 1996 ; 
|Gurevich et al. 1997] ). However, data arc too few to 
conclude that this interpretation is correct. We have 
examined spheroidal non singular isothermal density 
models, but it is still possible that the halo mass den- 
sity has a different r adial profile with a non fla t rota- 
tion curve (see, e.g., Wilkinson fc Evans 1999 ) which 
deserve some attention. 



The estimated ranges for a do overlap both for models 
with the same halo flattening q but different mass lim- 
its ()ii,n u ) and for models with the same mass limits 
but different halo flattening. This is encouraging since 
it means that our estimates do not depend on the halo 
flattening which is poorly constrained. Considering the 
whole set of models in Table 4, an estimate of a can 
be given in the range 0.27 -j- 1.90, which practically is 
the widest range of variability. Unfortunately it is not 
a narrow range because of the high errors connected 
to the low statistics. 

It is not very easy to compare this estimate with other 
results, since the present state of the studies on th e 
dark halo MF is still in progress ( Mera et al. 1998 ). 
Collecting together results on the nearby halo LF, 
based on geometric parallax determinations of high ve- 
locity faint stars in the solar neighbourhood, and those 
about the spheroid MF, Chabrier & Mera estimate 
that the halo MF is a power- law with a = 1.7 ± 0.2 



flChabrier fc Mera 1997| ) . We can see from Table 4 that 
the range for a obtained for the B2 models, that is for 
/i; = 0.01 and fi u — 1.0, are consistent with the value 
obtained by Chabrier & Mera (1997). This means to 
put a constrain on the low mass objects of the dark 
halo. Another important indication on the halo MF 
comes from studies of the stellar populations of the 
globular clusters. The MFs of different clusters ob- 
served with the HST have been determined recently 
by Chabrier & Mera (1997) who found that they are 
weakly dependent on the metallicity and may be mod- 
elled as power - law with slope a ~ 0.5 — 1.5 till 0.1 Mq. 
With this range for a also the A2 models are consis- 
tent. Note that we have also consistency with MACHO 
results. 

The estimated ranges for / are consistent 
with the recent estimate obtained by MACHO 



( |Alcock et al. 2000a|) using a maximum likelihood 
method and a different set of halo models. The best 
fit values fo are constrained in a very narrow range: 
from 0.19 to 0.24; the analysis of Table 4 shows that 
the total range for / is quite narrow as well: from 
0.10 to 0.35. This shows that it is not possible to 
increase the dark halo mass fraction flattening the 
dark halo. The same agreement with MACHO results 
is found with respect to the baryonic mass inside 50 
kpc (M^ r in Table 4. Our estimates range from 4.02 
to 15.12 x 10 10 Mq whilst MACHO estimates range 
from 6 to 13 x 10 10 Mq. Also the best fit values are 
in good agreement. These are encouraging results and 
lead us to be quite confident in our analysis. Once 
again, one should also consider different radial profiles 
for the mass density before drawing some definitive 
conclusion. 

Although not definitely conclusive, we obtained con- 
strains on the dark halo parameters on the basis of mi- 
crolensing observations for a wide class of spheroidal non 
singular isothermal models. We performed our analysis 
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on the basis of the microlensing data only in one di- 
rection, the LMC, and in fact it has been not possi- 
ble to determine uniquely the properties of the mod- 
els we considered. In order to investigate the shape of 
the halo and put more precise constraints on the rela- 
tive parameters, it will be necessary to have microlens- 
ing data relative to observations towards other direc- 
tions, for example towards globular clusters; they may 
be used as sources (Gyuk & Holder 1997) or as sites of 
lenses when observing towards SMC ( Jetzer et al. 1998| ). 
There are also other p ossible targets such as spiral arms 
( |Grcnachcr et al. 1999 ) which could be investigated. It 
would be also interesting to compare microlensing results 
on galactic models parameters with the ones obtained 
from the rotation curves of other spiral galaxies similar 
to our Milky Way, in particular those deduced from the 
universal rotation curve proposed by Persic, Salucci & Stel 
(1996) upon their analysis of a large homogenous sample 
of spiral galaxies. All that constitutes material for further 
work. 
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